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We provide an argument based on flux insertion to show that certain superconductors with a 
non-trivial topological invariant have protected zero modes in their vortex cores. This argument 
has the flavor of a two dimensional index theorem and applies to disordered systems as well. It 
also provides a new way of understanding the zero modes in the vortex cores of a spinless Px + ipy 
superconductor. Applying this approach to superconductors with and without time reversal and 
spin rotational symmetry, we predict the necessary and sufficient conditions for protected zero modes 
to exist in their vortices. 



Fermionic zero energy modes are predicted to exist in 
defects in a variety of condensed matter systems from one 
dimensional (Id) chains of polyacetylene to the vortices 
of certain two dimensional (2d) superconductors [2 [2]. 
These modes typically occur in systems with a bulk gap 
which are in a topologically non-trivial phase. The zero 
energy modes are often topologically protected in the 
sense that their existence is only tied to the particular 
topological phase of the system and to some global prop- 
erties of the system such as the winding number of the 
phase of the order parameter in a superconductor. 

Unpaired fermionic zero energy modes which are lo- 
calized in the vortices of certain superconductors such 
as a spinless px + ipy superconductor are of special in- 
terest. These modes, also called Majorana modes, have 
attracted a fair bit of recent theoretical attention due to 
the fact that they can obey non-abelian statistics and 
have potential applications in decoherence free quantum 
computation [3H6j . The existence of these modes has tra- 
ditionally been demonstrated through explicit analytic 
solutions of the Bogoliubov de-Gennes (BdG) equations 

[DiaiTHS]. 

In Id, Jackiw and Rebbi showed that the Dirac equa- 
tion interacting with a scalar field with a topologically 
interesting configuration such as a soliton has zero energy 
solutions whose number was related to an integer which 
characterizes the non-trivial topology [lOj . Similar zero 
energy states also occur in polyacetylene [5]. Recently 
many works have attempted to explain the Majorana 
modes in 2d by reducing the problem to a Id problem 

Since the vortex zero energy states occur in 2d sys- 
tems, a more general and inherently 2d explanation of 
the Majorana zero modes is desirable. In the current 
work, this is achieved through a flux insertion argument 
similar to Laughlin's gauge argument |14j . The argu- 
ment is not tied to any particular model and also applies 
to disordered systems where translational invariance is 
broken and is particularly useful in the context of recent 
efforts to find Majorana modes in systems other than a 
Px + iPy superconductor |15H18j . One of our main re- 
sults is a prediction about precisely what systems in all 
the different symmetry classes of the classification intro- 
duced by Altland and Zirnbauer[19l support zero energy 



states. 

Using a flux insertion argument, we first show that a 
certain class of insulators with tt flux inserted through a 
plaquette have exact zero energy eigenstates which are 
topologically protected in the sense that their existence 
is connected to the Hall conductance of the insulator. 
We then use these results to study superconductors with 
and without time reversal and spin rotational symmetry 
and predict precisely which systems have protected zero 
modes. 

The insulator Hamiltonians studied below have the 
same symmetries as the mean field Bogoliubov-de- 
Gcnnes Hamiltonians of superconductors. Though such 
insulators arc possibly of little physical significance by 
themselves, we study them since their spectra are identi- 
cal to that of the BdG Hamiltonians of superconductors. 
Further since they are insulators rather than supercon- 
ductors, they do not expel flux, making a flux insertion 
argument possible. Consider a 2d tightbinding insulator 
on a lattice with an even number of orbitals, 2s per site. 
The basis states are written in the form \i, a) where i is 
an index for the position and a for the orbital. We first 
study infinite sized systems where the energy of the topo- 
logically protected vortex states is exactly zero and later 
consider corrections that would occur in a finite sized 
system. 

The Hamiltonian can be written in the form Ji = 
J2i j ^i^ij'^j where Hij is a 2s x 2s dimensional ma- 
trix and 'iij ~ (-(/'j.q)"^, ^| = (V4/3) where i,j are the 
lattice position indices and a, ji the orbital and spin in- 
dices. The Fermi energy is set to zero and is assumed to 
lie in a gap. 

In the first instance, we study systems with neither 
spin rotational nor time reversal symmetry, though we 
will analyze systems with these symmetries later on. We 
further restrict our study to Hamiltonians H which posses 
a symmetry analogous to that of Bogoliubov de-Gennes 
Hamiltonians. In other words, we assume the existence 
of an anti-unitary operator, S*, such that 

sns-^ = -n. (1) 

Further, we assume that S acting on single particle po- 
sition eigenkets in the Hilbert space produces a linear 
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combination of kets at the same position, i.e., S^jS~^ = 
Uj'ij and Uj is a 2s x 2s dimensional unitary matrix. 
We may therefore conclude that the Hamiltonian of the 
system despite being that of an insulator, belongs to the 
symmetry class D in the classification introduced by Al- 
tland and Zirnbauer|19j. The symmetry under S also 
implies that U,HijU~^ = -H*j. 

Consider the effect of a vector potential which arises 
from flux insertion through an infinitesimal tube. The 
Hamiltonian of the system with a flux tube somewhere 
in the sample may be derived by multiplying each 
hopping term in the matrix Hij by the phase factors 

exp [i{e/h) dr. A] where the integral is along the hop- 
ping path and A is the vector potential. The Hamilto- 
nian of the system with the vector potential, which we 
call H{A) then has the following property: 

Sn{A)S-^ = 5(^*|i/,j*je*("/''^^'''''- (2) 

hi 

= -n{-A). 

Consider an infinite sample of the above system and 
let us adiabatically thread flux through an infinitesimal 
fiux tube through a plaquette at the center of the sample. 
If ^(0) is the Hamiltonian in the presence of a flux (/>, 
then on the basis of the above, we conclude that 

sn{(i))s-^ ^ -%{-(!)). (3) 

Thus if T-L{(t)) has an eigenstatc with eigenvalue E, then 
T-Li—cf)) has an eigenstate with eigenvalue —E. 

Now suppose the system has a quantized Hall conduc- 
tance of pe^ I^Trh. Then as a flux of 27rfi/e is adiabatically 
inserted through the flux tube, a total charge of pe flows 
in from infinity towards the flux tube |14j . In an infinite 
sample, the total spectral flow, i.e., the total number of 
states which cross the gap at the Fermi surface from be- 
low minus the number of states which cross it from above 
is equal to v [Till^U] . 

Let be the number of eigenstates of the Hamil- 
tonian whose energy is zero at the flux (j). A schematic 
example of the energy spectrum of extended states close 
to the Fermi energy as a function of flux is plotted in 
Fig. 1. The function n(0) is non-zero only for (j) e 
{0a, 0b, 0c, 0d} where (p^ = 2TT-(pa and 0^ = 2TT-(j>b. 

From Eq. 

n{2TT-(j))^n{(t)). (4) 

Apart from the states which traverse the gap and thus 
cross the Fermi energy an odd number of times, there 
might also be states which cross the Fermi energy an 
even number of times as shown in the flgure. 

Thus, 

J2 ?^(0z)=P + 2m, (5) 



where p is the total Chern number of the ground state. 
It follows from Eqs. Q and ^ that when p is odd, 

n(0) + n{2TTh/2e) ^ 2k + 1 

where k is some integer. 

At zero flux, the Fermi energy lies in a gap. Thus, 
n(0) = which in turn implies that n{2TTh/2e) is an odd 
integer when p is odd. Since the integer p is a topological 
invariant, which cannot change under small transforma- 
tions of the Hamiltonian, the zero mode is topologically 
protected. 

We now use this result to study superconductors. In 
the remainder of the paper, we are frequently going to 
consider Hamiltonians with 2 s degrees of freedom per 
site and we write such matrices in the form: 

/Mil Mi2\ 

Consider a tightbinding BCS Hamiltonian for fermions 
hopping on a lattice which can be written in the form: 

ns = E.,,*f^^.f*f where = (<^,„ is a 

Nambu spinor. Here, i, j stand for positions, and 7 is 
an index for the orbitals and spin which runs from 1 
to s. The matrix is the BdG Hamiltonian which 
is of the form of Eq. ^ with A/n — h, M12 — A and 
M22 — —h^ where h is the single particle Hamiltonian 
and A is the gap matrix. We can map the Hamiltonian 
Hs, to the Hamiltonian of an insulator, (which we call 
the associated insulator) given by Hi = J^ij ^l^Fj^ 

where '^j = {ipi^a, V'i,^)'^- 

We flrst study BdG Hamiltonians with neither time 
reversal nor spin rotational symmetry. Let us imagine 
inserting an infinitesimal flux tube containing a flux of 
2TTh/2e at the origin of our coordinates which is placed 
at the center of one of the plaquettes of the system. We 
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FIG. 1: A schematic plot of the energy, E versus (j), where 
(j) is the flux inserted through a plaquette of the insulator 
with Bogoliubov symmetry. Only states which lie close to 
the Fermi energy at = are shown. The dashed line is 
the single energy curve which traverses the gap as the flux is 
changed from to 2nh/e. 



3 



imagine that this flux does not leak into the rest of the 
superconductor and that the flux is inserted in such a 
way that the low energy configuration where the phase 
of the superconductor winds by 2tt around the flux tube 
is attained. 

In the presence of the flux, the single particle Hamil- 
tonian gets transformed as follows. Let r^, 9i be re- 
spectively the distance from the origin and the polar 
angle of site i in a coordinate system with the ori- 
gin located at the position of the flux tube. Then, 
/i^ = %e*(^/'*)/'*''-^ = /iye*(^--^^)/2, where A 

is chosen to be ||V(6'). The gap matrix transforms as: 
Aij A'^J = Ae^(^'+*j)/2 in the presence of the flux 
tube. 

The BdG eigenvalue equation in the presence of the 
vortex is thus: H'lp = Etp, where ip = {u, v)'^ and 
H' has the form of Eq. ^ with Mn = h',Mi2 = A' 
and M22 — Let Ui = Ui,Vi = e^^'Vi Then — 

(m, v)'^ satisfies the eigenvalue equation : H"ip — Eip, 
where H" can be written in the form of Eq. (pi) with 
Afn = h\Mi2 = A",M22 = -ihn,K'^ - /i,,e^W-«^)/2 
and A'/- = A,^-e*(^--«^)/2. 

We now replace the superconductor with the associ- 
ated insulator, "H/ with the same flux configuration, i.e., 
half a quantum of flux inserted at the origin, which lies at 
the center of a plaquette. It is easy to verify that T-Li{ti) 
can be written as '^\H^j'^ j with H" as given above. 
Further, T-Lj satisflcs the symmetry in Eq. ([T]) since it is 
derived from a BdG Hamiltonian and the analysis for the 
zero modes for insulators made previously can therefore 
be used. 

The necessary and sufficient condition for the existence 
of an exact zero-energy mode which is localized around 
the flux tube is therefore that the Hall conductance of 
Hi is pe^ /2Trh where p is odd. If this condition is met, it 
follows that there is a zero energy mode localized around 
the plaquette containing the tube. Since the Hall con- 
ductance is a robust topological invariant, the existence 
of the zero mode for the superconductor is also topolog- 
ically protected. 

A more realistic description of a vortex would include 
a finite region larger than a single plaquette where the 
magnetic field is non-zero rather than the situation con- 
sidered above where the flux is confined to a single pla- 
quette. In the general case, the gap parameter and 
the single particle Hamiltonians may be modeled as: 
A^ = A,,/((r, + r,)/2)e*(«'+^^)/2;/i^ = /i,,e'^-. Here, 
/(r) is some function which is zero at the origin and goes 
to one for distances larger than the coherence length and 
= / {e/h)dr.A where A is the vector potential which 
corresponds to a smeared out flux and which approaches 
the value {h/e)V{0)/2 for distances much larger than the 
penetration depth. 

The corresponding BdG Hamiltonian may therefore be 
written as Ji = Ho + V where Hq is the Hamiltonian 
which corresponds to the idealized vortex discussed above 
and 1^ is a local perturbation in the sense that Vij = 



sufficiently far from the vortex. If there is an odd number 
of Majorana zero energy states, then the Hilbert space 
of the idealized vortex has one unpaired fermionic mode, 
while every other eigenstate of energy, E, can be paired 
with a state of energy, —E. Any local perturbation can- 
not alter the Hilbert space structure. It follows that the 
perturbation, V, does not alter the existence of a zero en- 
ergy mode and that the zero mode persists in the more 
realistic configuration. 

We now intend to study superconductors with time re- 
versal symmetry using the same technique and therefore 
first study the corresponding insulators. Insulators with 
time reversal symmetry are classified by a Z2 invariant. 
Any Hamiltonian, H, of this class can be continuously 
deformed, without closing the gap or breaking time re- 
versal symmetry, to a Hamiltonian which is diagonal 
in the basis of the z-component of spin and can thus be 
written in the form H' — T-L'^ + H'^. The Hamiltonians 
H'^ and H'^ can be analyzed separately for studying their 
zero modes. The sum of the Hall conductance of H'^ and 
H'j^ must be zero. If laxy^nh/e^l of H'^ is an odd number, 
then the insulator belongs to the non-trivial class. By 
the previous analysis, in this case, when a flux of tt is 
inserted through the central plaquette the system has an 
odd number of Kramers pairs at zero energy, one mem- 
ber of each pair associated with "H^ and the other with 
"H^ . Kramers theorem then prevents a gap from opening 
up at TT flux when the Hamiltonian is continuously de- 
formed without breaking time reversal symmetry. This 
implies that the original Hamiltonian, H, must also have 
an odd number of pairs of zero modes at tt flux. When 
(Jxy{T-L'^) is an even multiple of e^/27r/i on the other hand, 
the number of pairs of zero modes must be even and is 
therefore not protected. 

Superconductors with time reversal symmetry fall in 
the class Dili of the classification scheme [T^. These 
superconductors are classified by a Z2 invariant [H] and 
can be mapped onto insulators with time reversal symme- 
try which satisfy Eq.([T]) and which have the same spec- 
trum, exactly as in the case of superconductors with- 
out time reversal or spin rotational symmetry. We have 
shown in the previous paragraph that when fiux is in- 
serted through a plaquette in these associated insulators, 
they have an odd number of Kramers pairs of zero en- 
ergy modes if and only if they belong to the non-trivial 
topological class. The zero eigenmodes of the equation 
that determines the spectra of the insulator with half a 
quantum of flux inserted through a plaquette may be re- 
lated to the zero eigenmodes of the superconductor with 
a vortex following the discussion after Eq. Thus, 
superconductors with time reversal symmetry have an 
odd number of Kramers pairs of zero energy modes in 
their vortex cores whenever the superconductors are in 
the non-trivial topological class. A pair of zero-energy 
modes may be regarded as a single Dirac mode. 

Hamiltonians of superconductors with spin-rotational 
symmetry but not time reversal symmetry fall in the 
class C of the Altland-Zirnbauer symmetry classes 
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[TO] , These Hamiltonians may be regarded as the 

sum of = (V-j, 'tlJi)H^{ip^, ipl)'^ and Hi = 

{tp^^, tlj^)Hi{tlji, '0j)^P^, where and have the 
form of Eq. Q with Afn = h, M12 = A, A/22 = ~h'^ 
and Mil = h, A/12 = —A, A/22 = —fi^ respectively. 

TABLE I: Conditions for superconductors in the various sym- 
metry classes to support protected zero modes in vortex cores, 
expressed as conditions on Hi, the insulator associated with 
the superconductor. The last column indicates whether there 
is a single protected Majorana mode(M) or a protected pair 
of modes(D). No protected modes exist for the class CI. 



Class 


Time-rev 


Spin-rot 


Condn. on Hi 


Mode 


D 


No 


No 


axy2Tvh/e^ — 2k — 1 


M 


C 


No 


Yes 


a^y2nh/e'^ = 2(2fc - I) 


D 


Dili 


Yes 


No 


non-trivial Z2 


D 


CI 


Yes 


Yes 







The spectra of and are identical. If (m, u)-^ 
written in the particle hole basis, is a zero mode of H^, 
then {u, —v)'^ is a zero mode of H^. The superconductor 
may be mapped onto an insulator, "H/, which is the sum 
of two single particle Hamiltonians, Hi^^ and "H/,^, and 
which may be regarded as separate systems. The con- 
dition that in the presence of a vortex, the matrix 
has an odd number of zero modes is, as deduced in the 
study of Hamiltonians of class D, that the Hall conduc- 
tance of the corresponding insulator, axy{'Hi^'f-)2Trh/e^ is 
an odd integer. The net Hall conductance, cFxyiHj) is 
twice that of ///,^ and is thus always an even integer. 
Thus, when the Hall conductance of Hi has the form 
2pe^/2Trh where p is an odd integer, the system has an 
odd number of pairs of zero modes in its vortices, while 
when p is even, the system has an even (or zero) number 
of pairs of zero modes. 

Superconductors with both time reversal and spin ro- 
tational symmetry, which belong to the class CI, may be 
regarded as belonging to the trivial Z2 class of supercon- 
ductors with time reversal symmetry. These supercon- 
ductors thus have no topologically protected zero modes 
in their vortex cores. 



Our results are summarized in Table H] In the notation 
used in Ref. [22], superconductors in the class D have 
protected Majorana modes in their vortex cores when the 
integer invariant for this class is odd, superconductors in 
class C and DIII have protected Dirac modes when the 
invariants for these classes are odd and non-trivial respec- 
tively. Most continuum models can be simulated to an 
arbitrary degree of accuracy by a series of lattice models. 
Since the results derived above are not limited to a par- 
ticular tightbinding model, one expects that the analysis 
presented above extends also to continuum models 24 . 

Our analysis thus far has been restricted to a single 
vortex in an infinite superconductor. A finite supercon- 
ductor with a vortex is topologically equivalent to a sys- 
tem with two vortices and periodic boundary conditions. 
One may then start with the associated insulator with 
no flux inserted and gradually insert flux into two tubes 
in such a way that the flux entering one tube is equal to 
the flux leaving the other leaving the total flux through 
the system zero. At 27r/i/2e flux, the states localized at 
the two flux tubes will in general hybridize, giving rise to 
a finite splitting. The magnitude of the splitting when 
the two vortices are a certain distance d is proportional 
to the overlap between two zero-energy eigenstates of the 
infinite system placed at the same distance and therefore 
falls exponentially as the distance between the vortices. 

Finally, we note that the arguments made above would 
are also applicable in the case when there is a mobility 
gap in the absence of flux rather than a gap to all states 
since the states within the mobility gap have zero Chern 
number and their energies return to their initial values 
when the flux inserted varies from to 27r/i/e. 

In summary, we have provided a simple and general 
argument which shows that certain topological classes 
of superconductors have topologically protected, robust 
zero modes, which can either be unpaired Majorana 
modes, or come in pairs. We applied this analysis to 
the various symmetry classes of superconducting Hamil- 
tonians. 
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